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SECONDARY INSTABILITY OF COMPRESSIBLE BOUNDARY LAYER 


TO SUBHARMONIC THREE-DIMENSIONAL DISTURBANCES 
Nabil M. El-Hady 

Department of Mechanical Engineering and Mechanics 
Old Dominion University, Norfolk, Virginia 23508 


ABSTRACT 


Three-Dimensional linear secondary instability theory is extended for 
compressible boundary layers on a flat plate in the presence of finite amplitude 
Tollmien-Schlichting (T-S) wave. The focus is on principal parametric resonance 
responsible for strong growth of subharmonics in low disturbance environment. 


I. INTRODUCTION 


Whatever theoretically or experimentally, the compressible stability theory 
lack a firm connection with boundary-layer transition. "There is little doubt 
that transition is preceded by linear instability in many instances, but the way 
these individual unstable waves act, alone or in combination , to trigger the 
transition process is not known" [1]. The nonlinear theories and the secondary 
instability that are much a prominent feature of incompressible stability theory 
do not exist for compressible boundary layer. 


1 


Recently, Erlebacher and Hussaini [2,3] by using direct simulation of 
parallel boundary layer, generated numerically a high Mach number vortical 
structure (peak-vally) similar to that observed and computed for incompressible 
flows. 

Stimulated by this work, we formulated the secondary three-dimensional 
instability problem for compressible boundary layers. In this paper, we 
investigate theoretically the effect of finite amplitude two-dimensional (2D) 
wave on the growth of three-dimensional (3D) perturbations in compressible 
boundary layers. Hence, this paper covers only a range of Mach numbers up to 
the transonic, where the critical primary disturbance is 2D. For supersonic 
boundary layers, the critical primary disturbance is 3D, or a 2D second mode. 
This range of Mach numbers is currently under investigation. An analysis 
similar to that of Herbert [4], Herbert et al [5] and Nayfeh [6] is followed 
here but spatial stability is considered for both primary and secondary 
instabilities. 

The primary instability leads to the growth of T-S waves and a streamwise 
almost periodic modulation of the flow. We study the linear stability of this 
flow with respect to spanwise periodic 3D disturbance. Floquet theory gives as 
a solution to the stability equations, all various types of resonance. We 
consider the case of principal parametric resonance responsible for strong 
growth of subharmonics in low disturbance environment. 


II. PRIMARY INSTABILITY AND THE BASIC FLOW 
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We cosider the primary instability of the 2D compressible boundary layer on 
a flat plate, with respect to 2D quasi-parallel spatially growing T-S 
disturbances. The critical disturbance is 2D up to the transonic range of Mach 
numbers [7]. The flow field is governed by the 2D compressible Navier-Stokes 
equations , the energy equation and the state equation. Dimensionless quantities 
are introduced by using the reference velocity U r and the referece length 
L = x/ Ur , so that Reynolds number is given by R = /u,* x/m , where x measures 
the distance from the leading edge of the plate, and M is the fluid kinematic 
viscosity. The thermodynamic and transport properties of air (treated as 
perfect gas) are made dimensionless using their corresponding freestream values. 

At sufficiently large distance from the leading edge, primary instability 
of the compressible laminar flow occurs with respect to T-S disturbances. These 
disturbances take the traveling wave form, 

(1) (up ,vp ,pp ,tp ,pp,pp ) = A 0 [ u(y),v(y),p(y),t(y),p(y),p(y) ] 

exp [i(Jo< dx -C6t)] + c.c. 

Where the y-axis is defined normal to the flat plate, u and v are the 
disturbance velocities, p is the disturbance pressure, t is the disturbance 
temperature, p is the disturbance viscosity, and p is the disturbance density. 

For the spatial stability analysis o< is a complex wavenumber given by 
o< =<Xr + and (0 is a real disturbance frequency, and c.c. denotes complex 
conjugate terms. The eigensolutions u,v,p and t are governed by a six-order 
system of equations that is given in reference (9). The density disturbance is 
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related to the temperature and pressure disturbances through the state equation, 


(2) V P = p T + p t 

while the viscosity disturbance is related to the temperature disturbance by 

(3) p - ( dp /dT ) t - p t 

where p and p are the mean-flow density and viscosity respectively,^ is the 

ratio of specific heats and M^is the freestream Mach number. The six-order 

system is numerically integrated as initial value problem using a combination of 

shooting [8] and Newton-Raphson iteration technique that employs a Gram-Schmidt 

orthonormalization procedure. The solution of this eigenvalue problem is 
calculated by the author elsewhere [7,9]. 

The linear stability theory of primary instability provides for a given tO 
and R. Then the integration of the growth rate — gives the amplification 
factor , 

( R 

(4) In ( A / A. ) - -2 J dR 

R-op 

where A 0 is an arbitrary initial amplitude of the primary instability at R op (R 

AS ' 

where the onset of the primary wave) . The eigensolutions may be normalized 

such that A measures directly the maximum r.m.s. value of the streamwise 

disturbance, that is 

2 

(5) max (u(y)| =1/2 

0*y<oo 

Since the primary instability of boundary-layer flows is induced by viscosity, 
the growth rates and amplification factors here are typically very small. 
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Now the basic state under consideration is composed of the sum of the 2D 
compressible mean flow and a finite amplitude A of the primary T-S wave, that is 

(6) (u^ ,v^ ,p fc ,t^ ) = (U,0,P,T) + A [ u(y),v(y),p(y),t(y) ] exp(i0) 

Where 

A = A o exp(- ^ Otfj dx) , assumed constant, and 


e 



dx - t 


III. SECONDARY SUBHARMONIC INSTABILITY 


We consider the 3D quasi-parallel spatial subharmonic instability of the 
basic state given by (6). The finite amplitude primary wave acts as a 
parametric exitation on the secondary instability. Following the analysis of 
Herbert [4,5] , Nayfeh [6], and El-Hady [10] we apply Floquet theory and express 
the secondary wave using the normal mode concept, 

A A A A /V \ 

( u s »v s ,p s ,t ft ) = exp(^Y dx) [u(y) , v(y) ,p(y) , t(y) ] exp(^ i0) cos pz 


(7) 


+c.c. 


w s = exp(^V dx)w(y) exp(4 i0) sin pz + c.c. 


A 

where p is a spanwise real wavenumber, and Y = i)f| is a characteristic 
exponent. The spatial growth rate of the secondary wave is given by Yf , while 
Y; can be interpreted as a shi&t in the streamwise wavenumber. In our 
calculations, we consider only the case of Y{ = 0, that is the secondary wave is 
perfectly synchronized with the basic state. 
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The secondary wave (7) is superposed on the basic state (6) and the result 
is substituted into the dimensionless Navier-Stok.es equations. The meam flow 
plus the 2D T-S quantities are substracted, and the resulting equations are 
linearized in the secondary disturbances. Then one obtains an eigenvalue 
problem that can be written as 

-I si 2 " ' 

(8) T ( g Z, - T DT Z 3 + DZ + p Z^ ) + G ( J M*, Z g - T Z-j ) 

+ A (aj ) = 0 

(9) { G - p R* [(nri-2) i ~ p 2 ] } Z v - R ' D( f Zj ) 

-l - 

+( T DU - Dyi R g ) Z 3 - (nri-1) ji R g (DZ^ + p Z^ ) 

•I ^ a —I 2 

+ g Z 6 - R D(p DU) Z-j - yi R DU Z ft + A (a^) + A (b 2 ) = 0 

(10) [ G - p r"‘ ( £ + f ) ] Z 3 - (nri-1) p r'( g Z 2 + p Z s ) 

•I * -l — 

-(m+2) R D(p DZ 3 ) - m R Dp ( g Z , + p Z^ ) + DZ g 

A -I Z 

- p R DU g Z-j + A (a 3 ) + A (b^) = 0 

- -I - -I „ 

(11) (nri-1) p R p ( g Z, + DZ 3 ) + Dp R p Z3 

[ G - p R*‘ g 2 + (nri-2) p R l (3 Z ] Z 4 - R ' D( p Z g ) 

Z 

- P Z 6 + A (a 4 ) + A (b^) = 0 
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( 12 ) 


* - -\ A 

-2M p R DU ( Z 2 + g Z 3 )+ T DT Zj - M T G Z 4 

- a -l 

+{ G - R Pr [p (g - p ) + D(p DT)] - M R (DU) p } Z ^ 

—I —l _ — Z 

- R Pr ( 2 Z ft Dp - p DZ g ) + A (a s ) + A (b^) = 0 

where 


Z. 

A 

= u , 

A 

= Du , 

A 

Z s = v 

> 

A 

Z 4 = w 

z s 

A 

- Dw , 

A 

z 6 “ P > 

A 

z 1 = t 

> 

A 

Z g = Dt 


g = ( + £ i w r ) 

G = t'' ( gU - £ it*>) 

A v/ 2- 

M = (Y- 1 ) M w 


also Pr is Prandtl number, D=d/dy, and m=2(e-l)/3 is the ratio of 
coefficients of viscosities, where e=0 corresponds to the Stokes hypothesis. 
The boundary conditions are 


(13) 


Z | = Zj = Z^ = Zj = 0 at y= 0 


Z | , Z j , Z^ , Z i -» 0 as y ■>« 


the 
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Equations (8)-(12) govern the secondary 3D instability of compressible 2D 

flows. They represent the mass, x-momentum, y-momentum, z-momentum, and energy 

2 

equations respectively. The coefficients of A and A ( a; , bj , i=l-5 ) are given 

2 

in Appendix A. The coefficients of A in the mass equation and those of A in 
all other equations are new in this system of equations in contrast to the 
secondary instability equations for incompressible flows. These coefficients 
are present mainly due to the density disturbance of both primary and secondary 
instabilities. In the incompressible limit eqations (8)— (12) reduce exactly to 
that given by El-Hady [10] • When A=0,the system of equations (8)— (13) govern a 
primary subhamonic 3D wave. 

The system of equations (8)-(12) can be written as eight complex equations 
in the form, 

DZ t - Z z = 0 
DZ^ + c 0 DZ^ + Cj - 0 

DZj + c^ DZ^ + C£ = 0 

DZ /% " Z r = 0 

(14) * 5 

DZg + DZ«g "f Cg ~ 0 

DZ^ + c^ DZ^ + c^ = 0 

DZ ^ - Zg = 0 

DZ r + c g + c <\ = 0 

Where the c's are quadratures in the primary and secondary disturbance 
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quantities, and the overbar indicates a complex conjugate. Equations (14) are 
numerically integrated as initial value problem from y ® (edge of the 
boundary layer) to the wall. The eigenvalue search used a Newton-Raphson 
iteration technique to satisfy the last boundary condition at the wall. A well 
tested code SUPORT [8] is used which is coupled with an orthonormalization test 
based on the modefied Gram-Schmidt procedure to overcome the stiffness of the 
integrated system of equations. 

The linear stability theory of the secondary instability provides Y for a 
given /3 and R. Then the integration of the growth rate Yr gives the 
amplification factor, 

r R . 

(15) In ( S / S 0 ) = 2 ) } r dR 

R oS 

where S o is an arbitrary initial amplitude of the secondary instability at R oS (R 
where the onset of the secondary wave) . The secondary subharmonic instability 
is believed to originate from a strong mechanism of combined tilting and 
stretching of the vortices such as the case of incompressible flows [10], thus 
the growth rates and amplification factors are expected to be large as they 
occur on a convective length scale. 


IV. RESULTS AND DISCUSSION 


For the incompressible limit (M^ ■ 0), our compressible secondary 
instability code gives results that are in full agreement with those obtained by 
Herbert et al.[5] and by El-Hady [10]. All results reported here are for the 
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nondimensional frequency F * 6d / R * 60E-6 , that remains fixed as a wave of 
fixed physical frequency travels downstream. 

At 0, a primary 2D instability grows between R 0 p “ 554 and Rjp = 1052 
(first and second neutral points) reaching a maximum amplification factor of 
A/AO *41.68. As Mach number increases, the growth rates of the primary 2D 
waves decrease as shown in Fig.(l). Also the first and second neutral points as 
well as the streamwise location where the maximum growth rate occurs shift to 
the left, they occur earlier upstream as Mach number increases. Fig. (2) shows a 
reduction in the amplification factor of the primary 2D instability as Mach 
number increases. 

Almost in the same region where the primary 2D waves are growing, a broad 
band of spanwise wavenumbers of primary 3D subharmonic waves are subject to 
amplification. Fig. (3) shows the growth rate curves of these primary 3D 
subharmonics (F * 30E-6) for different Mach numbers. These curves possess the 
same features of the primary 2D waves of having lower growth rates that shift to 
lower R as Mach number increases. Fig. (4) shows the amplification factors of 
these primary subharmonics. Both Figs. (3) and (4) suggest that the growth rates 
and amplification factors of the primary subharmonics are typically so small to 
bear any resemblence to experimentally observed transition. 

However, the growth rates shown in Fig. (5) as function of the spanwise 
wavenumber B =10E+3 p/R for secondary 3D subharmonics, are much larger than 
those for primary 2D waves or primary 3D subharmonics. This strong growth is 
due to the parametrical exitation by the finite amplitude primary 2D wave. 
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Fig. (5) is calculated at R=850 and a primary 2D amplitude A =.01, for different 
Mach numbers. For high spanwise wavenumbers (B^.225), compressibility appears 
to have a stabilizing effect on the secondary subharmonic instability. In other 
words, the secondary subharmonic instability is largest at 0, and decreases 
as Mach number increases. But for low spanwise wavenumbers (B<.225) the effect 
of compressibility depends on the value of the spanwise wavenumber. 

At fixed F and R, results for different Mach numbers show a destabilizing 
effect ( higher growth rates for the secondary subharmonic ) as the primary 2D 
amplitude A increases. At R=1050, Fig. (6) illustrates the destabilizing effect 
of increasing A for M^O.8 compared with that for M^O. Fig. (6) reveals that 
the influence of compressibility on secondary subharmonics is not only function 
of the spanwise wavenumber, but also is function of the primary 2D amplitude. 
While at A=.01 compressibility is stabilizing for large values of B and 
destabilizing for small values of B, it has a destabilizing influence on the 
secondary subharmonics at A=.002 for all range of spanwise wavenumbers. The 
figure also shows that at very small amplitudes, considerable growth rates exist 
in a small band of wavenumbers, that extends to larger values as the amplitude 
increases. 

Fig. (7) shows the effect of Reynolds number R on the growth of the 
secondary subharmonics at M^O.8. As R increases , an increase in the growth 
rates exist at fixed F and A. 
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To evaluate the overall effect of compressibility on the secondary 
subharmonics, we can combine the influence of increasing the amplitude A of the 
primary 2D wave and increasing R for various Mach numbers at fixed F. For 
comparison purposes, the amplification factor of the subharmonic is calculated 
using equation (15) from R oS ( onset of the secondary subharmonic instability) to 
any R downstream. For different Mach numbers, Fig. (8) shows the variation of 
the growth rates of the secondary subharmonics with R at a spanwise wavenumber 
B=.15. The initial primary 2D amplitude used for these calculations is A=.001. 
Fig. (9) shows the amplification factors decreasing from about 29 at M^O. tO 12 

at M =0.8. 

00 

Figs. (8) and (9) indicate that while the growth rates and amplification 
factors of the secondary subharmonics are decreased by increasing Mach 
number, the onset of the secondary instability is almost not affected. This is 
probably due to the combined effect of the upstream shift of the onset of the 
primary 2D instability, as well as the growing of the primary amplitude. 


V. CONCLUSION 


We formulated the secondary three-dimensional instability problem for 
compressible boundary layers. The effect of finite amplitude two-dimensional 
T-S wave on the growth of three-dimensional subharmonics is investigated for a 
range of Mach numbers from 0 to 1.2. Numerical results for F =60E-6 show that 
the local (at fixed R) effect of compressibility on the secondary subharmonics 
may be stabilizing or destabilizing depending upon their spanwise wavenumbers. 
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as well as the finite amplitude of the primary 2D wave . However , the overall 
effect of increasing Mach number is a reduction in the growth rates and 
amplification factors of the secondary subharmonics, almost with no change in 
the streamwise location where this instability sets in. 
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APPENDIX A 

a, = f(g Z, +p ) + D(f Z 3 ) +g«F+ D(v F) 

by = 0 

_-l a _ _ a Z _ 

a£ “ [ G f + T gu-m 2 Rg(g + i °< ) t + R p t ] %\ 

—\ — | A — -l —l A — 

+ [ T V - R D( ji t)] Zj + [ f DU + T Du - R g D( |i t)] Zj 

A— A — A _ — 

- R t DZ 2 - R ( mg + g ) t DZ a + R p ( m | g + im«) t Z,q 

a _ _ 2 - 

- R [m 2 (X (i g +ot ) u + (im > 0<+ m g) Dv + D u ] Z*j 

- R ( Du + i« v ) D( p Z ? ) + [ i(« U -40 ) + v DU ] F 

=> [ (g + ) fu + gfu] Z i + ( fv + fv ) Z a + ( fDu + fDu ) Z 3 

+ ( v Du + v Du ) F 

— —I — A _ 2 ^ — 

a g = [ f G + T (gu + Dv )-R( g +i g«-p )t]Z 3 

+ [ t“‘ v - m z R** D( p t) ] DZ 3 - m 2 R t D 2 Z 3 + i( «U - ^ ) v F 

+ [ iT~ *<X v -mR 1 D( p t) i ] Z, - R (a, g + i»t) t Z 2 - mR p D( p t) 
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IS I 


A _ A _ Z 

Rp t Z^ - { R [ (im®<+g) Du+ (i g -<X) v + m^Dv ] 

”1 A - A — 

-R Dp ( im« u + mjDv )} Z-^-R ( it « u + t 2 Dv ] 


b 3 = io<( f v-f v ) Z, + [( f u + fu) g + fDv + fDv ] Z^ 


+ (fv + fv) DZ 3 + [ 2vDv + iot( vu^-vu )] F 


— -t — A _Z _ 

a^=[fG + T g u - R ( g + ifi* g - m p ) t ] Z^ 

—l -I A — A _ — — 

+ [T v-R D( p t) ] Z 5 - R [ t DZj - p t ( g + m g ) Z, 

- m, p t DZ 3 - m p ( i <X u + Dv ) Z^ ] + R - ' p D( f t) Z 3 


’ ( f U + f u ) g + ( f T + f T ) Zj 


, -- -I - AA -| A _ 2 _ 

a g - {£ G + T g u - 2M R DU (Du + i et v) - Pr [R t ( g + 2io<g - 
+ R* D( y Dt )]} Z 7 + [ i («u -tO) t + vDT ]F-MguZ 6 


A -T -l -I A A -** A _ 

- M v DZ 6 + [ T v - Pr ( 2 R Dt + t R Dp ) ] Z g - R Pr t DZ 

-> A A _ 

+ [ i °< ( T t - M p ) - 2 H p R g ( i « mg u - m Dv ) ] Zj 


A “I A — „ — 

- 2 M R [ p (i« v + Du)+ p t DU ] ( Z £ + g ) 

—I A 

+ ( T Dt + f DT - M Dp ) Z 3 


A —l _ _ 

-2MR ptCimOtu + m^Dv) DZ^ + m p ( iocu + Dv)Z^ ] 


2 2 
* - P ) 
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— — A A 2 _ __ _ 2. _ 

b_={ (fu + fu)g-2MR[o< vv + iot(vDu-vDu)+m 2 o<uu 


+ m ot ( uv + uv ) + m £ Dv Dv + Du Du ] Z. 


+ [io<(ut-ut) + vDt + vDt]F + (fv + fv)Z. 


+ [ic*(tf-tf)-2MR im 2 e<g(tu-t'u)]Z, 

A A _ _ 

- 2 M R ( t Du + t Du + i« ( tv-tv) (Zj+gZj) 


+ ( f Dt + f Dt ) Z* - 2M R i- ( t Dv + t Dv ) DZ; 


-I 

f-T ( 0 p “ T t) 


-1/2. -I 

F = T (J M^Z^ - T Z, ) 


* A -1 
R = p R 


m, ■ ra+1, = nri-2 


and the overbar denotes complex conjugate. The overbar in p and the subscript 
r in are omitted for simplicity ( p is p and Ot is Off). 
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Fig.l Growth rates of the primary 2D wave (F=60E-6) at 
different Mach numbers. 
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Fig. 2 Amplification factors of the primary 2D wave 
(F=60E-6) at different Mach numbers. 
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Fig. 3 Growth rates of the primary 3D subharmonic 
(F=30E-6 and B=.12) at different Mach numbers. 
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Fig. 4 Amplification factors of the primary 3D 
subharmonic wave (F=30E-6 and B=.12) at different Mach 
numbers. 
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Fig, 5 Growth rates of the secondary 3 
R=850 and A=.01 as function of the sp 
for different Mach numbers. 
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Fig. 6 Effect of the primary amplitude A on the growth 
rates of the secondary 3D subharmonic at R=1050, for 
Mach numbers 0 and 0.8. 
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Fig. 9 Variation of the amplification factors with R of 
a secondary 3D subharmonic at B=.15 for different Mach 
numbers . 
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